{

 TFile *f = new TFile("qcd_mm_d0sig_5.root");
 TH1D *hqcd_data = f->Get("data10_njets_qcd_estat_j0"); hqcd_data->SetDirectory(0);
 TH1D *hqcd_mc   = f->Get("mc_qcd_njets_j0");         hqcd_mc->SetDirectory(0);
 TH1D *hqcdbb    = f->Get("mc_qcdbb_njets_j0");       hqcdbb->SetDirectory(0);

 hqcd_mc->Add(hqcdbb);

 Double_t qcd_data_est[5] = {0,0,0,0,0};
 Double_t qcd_mc_est[5]   = {0,0,0,0,0};
 Double_t qcd_est_err[5]  = {0,0,0,0,0};

 qcd_est_by_jets(hqcd_mc,   qcd_mc_est,   qcd_est_err); 
 qcd_est_by_jets(hqcd_data, qcd_data_est, qcd_est_err);

 for (int i=0; i<5; i++) {
     Double_t r = qcd_data_est[i] / qcd_mc_est[i];
     cout << i << " " << qcd_data_est[i] << " " << qcd_mc_est[i] << " " << r << endl;
 }

}
